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1. INTRODUCTION 



A smoking-gun signature of the dark matter (DM) signal in DM direct detection exper- 
iments is the annual modulation of the event rate. This arises because the earth rotates 
around the sun, while at the same time the sun moves relative to the DM halo [1, 2]. 
Currently, two experiments report annual modulation of the signals, DAMA with 8.9a sig- 
nificance [3], and CoGeNT with a significance of 2.8a [4]. An important question is whether 
the observed modulation rates are consistent with the interpretation that the signal is due to 
DM. Ideally, one would like to answer this question without requiring a detailed knowledge 
of the DM halo. As we will show below, we come very close to this ideal. 

The DM scattering rate in the detector is determined by the particle physics properties of 
DM and by the astrophysical properties of the DM halo. The latter carry large uncertainties. 
When presenting results of DM direct detection experiments, it has become customary to 
use a "standard" Maxwellian DM velocity distribution in order to show the constraints on 
DM mass and scattering cross section. This is very likely an oversimplification, with N-body 
simulations indicating a more complicated structure of the DM halo, see e.g., [5]. It is well 
known that the annual modulation signal is sensitive to such deviations from the Maxwellian 
halo [6-8]. DM streams are particularly relevant for the annual modulation signal and can 
lead to large effects. These can, for instance, be comparable in size to the modulation effect 
due to the whole Maxwellian halo [9-12]. 

In the following we will derive constraints on the amplitude of the annual modulation 
using measured unmodulated rates. The derivation of these bounds is made possible by 
the fact that the velocity of the earth rotating around the sun, Vg ~ 30 km/s, is much 
smaller than the typical velocity of DM particles in the halo, (v) ~ 200 km/s, and we can 
expand the velocity integral to first order in Ve. The derived constraints on the annual 
modulation then serve as a consistency check for the hypothesis that the annual modulation 
is due to the DM signal. The check that we provide is (almost) independent of DM halo 
properties. Previously halo independent methods to interpret DM direct detection data 
have been developed and applied in [13-18], with different goals than ours. In ref. [14] 
a method for determining the DM mass independently of the halo was developed, while 
refs. [13, 18] focused on extracting halo properties from the data. Compatibility studies 
of different experiments without adopting assumptions on the halo have been performed 
in [15-17]. Our goal is to obtain a consistency check between the modulation amplitude and 
the unmodulated signal within a given experiment. 

The structure of the paper is as follows: in sec. 2 we set up the notation, followed in 
sec. 3 by a derivation of the general bound on the modulation signal under the very mild 
assumption that the properties of the DM halo in our vicinity do not change on time scales 
of months. In sec. 4 we impose further symmetry requirements on the DM velocity distri- 
bution and obtain successively more stringent bounds. In sec. 5 we then demonstrate the 
usefulness of our bounds by applying them to the modulation signals reported in CoGeNT 
and DAMA within elastic spin-independent scattering. We find that our method leads to 
non-trivial restrictions on DM interpretations of the modulation signal in CoGeNT whereas 
the DAMA signal satisfies our bounds for the relevant DM masses. We summarize in sec. 6. 
In appendix A we consider the Maxwellian DM halo and show that the expansion in Ve is 
rather accurate in this case. Furthermore, we illustrate the various bounds compared to the 
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modulation signal expected for the Maxwellian halo. Appendix B provides technical details 
necessary to translate the bounds derived for the halo velocity integral into observable event 
rates. Supplementary material related to the treatment of an unknown background is given 
in appendix C. 



2. NOTATION 

The differential rate in events/keV/kg/day for DM x to scatter off a nucleus (A, Z) and 
depositing the nuclear recoil energy Enr in the detector is 

R{E^r.t) = [ d\^vU.,{.^t). (1) 

Here ~ 0.3 GeV/cm^ is the local DM density, rriA and are the nucleus and DM masses, 
aA the DM-nucleus scattering cross section^ and v the 3-vector relative velocity between 
DM and the nucleus, while v = |v|. For a DM particle to deposit recoil energy Enr in the 
detector a minimal velocity Vm is required, restricting the integral over velocities in (1). For 
elastic scattering it is 



lmAE„ 



2n^ • 

where p^-A is the reduced mass of the DM-nucleus system. 

The function /det(v, t) describes the distribution of DM particle velocities in the detector 
rest frame with /det(v, t) >0 and / (Pvfdet{y, t) = 1. It is related to the velocity distribution 
in the rest frame of the sun by 

/det(v, t) = /sun(v + Ve(t)) , (3) 

where Ve(t) is the velocity vector of the earth, which we write as [9] 

Ve(t) = fe[ei sin X{t) — 62 cos X{t)] (4) 

with Ve = 29.8 km/s, and X{t) = 27r(t — 0.218) with t in units of 1 year and t = at 
January 1st, while ei = (-0.0670,0.4927,-0.8676) and 63 = (-0.9931,-0.1170,0.01032) 
are orthogonal unit vectors spanning the plane of the earth's orbit, assumed to be circular. 
Similarly, the DM velocity distribution in the galactic frame is connected to the one in the 
rest frame of the sun by /sun(v) = /gai(v + Vgun), with Vgun ~ (0, 220, 0) km/s + Vpec and 
Vpec ~ (10,13,7) km/s the peculiar velocity of the sun. We are using galactic coordinates 
where x points towards the galactic center, y in the direction of the galactic rotation, and z 
towards the galactic north, perpendicular to the disc. As shown in [7], eq. (4) provides an 
excellent approximation to describe the annual modulation signal. 



^ Throughout this work we assume that DM is dominated by a single particle species. A generalization to 
mult i- component DM is straightforward and will be pursued in future work. 
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In the following we consider the typical situation, where the differential cross section is 
given by 

da A tua Tp2 



dEnr '^IJ'Ia'"^ 



a\F\Enr) , (5) 



where a\ is the total DM-nucleus scattering cross section at zero momentum transfer, and 
F(Enr) is a form factor. The event rate is then given by 

R{Enr,t)^CF\Enr)v{Vm:t) with C - (6) 

and the halo integral 

VM^I dhl^^. (7) 

Here and in the following Vm and E^r are related by eq. (2). This formalism covers a wide 
range of possible DM-nucleus interaction models, including the standard spin-independent 
and spin-dependent scattering. The results we derive below apply to all the cases where 
daA/dEnr oc in which case the halo integral (7) is obtained. The arguments can easily 
be generalized to a non-trivial dependence of the interaction, which would introduce an 
additional Enr dependent function in (6) but would not change (7). Our results do not 
apply, however, to a non-standard velocity dependence of the cross section, which would 
modify eq. (7). 

The rate will have a time independent component and an annually modulated component 
which we define as 

t) = R{Enr) + SR{E,,r, t) = C F\Er,r) [f]{Vm) + H^m, ^)] • (§) 

Below we will be specifically interested in purely sinusoidal time dependence with period of 
one year, in which case we can write 

SR{Enr, t) = AR{Enr) COs27r[t - to(£'nr)] , /gx 
SviVm, t) = Ar,{Vm) COS 2'!r[t - io(-E^nr-)] ■ 

The peak of the annual modulation occurs at to, which in general depends on Enr (or 
equivalently on Vm)- The modulation amplitudes of the event rate, A^, and of the halo 
integral, Ajj, are related through Aji{Enr) = CF'^{Enr)Ari{vm) > 0. We will first derive 
bounds on A^ in terms of the time averaged value of the halo integral fj. In section 5 we will 
then translate the bounds into constraints involving the observable quantities Ar and R. 



3. GENERAL BOUND ON THE ANNUAL MODULATION AMPLITUDE 

Assumption 1: We assume that the only time dependence comes from Ve{t) and there 
is no explicit time dependence in /sun- 

This assumption implies that the halo is spatially constant at the scale of the sun-earth 
distance and also constant in time on scale of months. Under this assumption we can derive 
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a general bound on the annual modulation by expanding eq. (7) in Ve/vm <^ 1 ^- Using for 
short / = /sun, we have 

Vivm,t)=[ d'^J^i (10) 

73„./(v) , /"^3„,^^,,aV- V^(t) ^^^^^ ^ ^_ ^ ^ ,nf..2/..2 

,3 



V>Vm 



(11) 



where the first term in (11) gives the time independent part of the DM scattering signal. In 
polar coordinates the time independent halo integral is then given by 

f]{vm)= [ d^v^^= I dvv I dip I dcos^f{v,'d,ip). (12) 

J V>Vm ^ J Vm Jo J —I 

The second term in (11) corresponds to the time dependent part of the DM scattering signal, 
with Ve(t) given in (4). Expanding to linear order in thus leads to an annual modulation 
signal that has a sinusoidal shape. One can check experimentally for the convergence of 
the expansion by searching for higher order terms, oc sin^[27r(i — to)], etc. Note that the 
expansion is in v^^/v, where v > so that the accuracy is typically better than 0{ve/v„i). 
In appendix A we demonstrate explicitly that this expansion is rather accurate in the case 
of a Maxwellian halo. 

By expanding the velocity integral in the small quantity we assume that /(v) is 
smooth on scales < v^- Hence, our bounds do not apply in situations where /(v) has strong 
structures at scales smaller than v^- An example would be a very cold stream with velocity 
^stream a^d a dispcrsiou smaller than fg. In such a case the expansion will not be accurate 
for Vm in the range \vm — fstreaml ^ feSina, though it may still work to good accuracy for 
Vm outside this range. 

The time dependent component in (11) has two contributions 

SviVm: t)^ j dhfiv)^^^^ [e{v - Vj - S{v - vjv^] . (13) 

The term with the 0-function comes from expanding the denominator in eq. (10) and involves 
a volume integral over the region v > Vm- The term with the 5-function comes from taking 
into account the effect of Ve on the integration boundary, and the 5-function indicates that 
the argument has to be evaluated at the surface v — Vm- Let us treat the two terms 
separately and define 

^g{Vm)g{Vm) = J d^V f (y) ^d{v - Vm) , (14) 

Mvm)G{vJ = J dhf{v)^e{v - vj . (15) 



^ Typically is sensitive to both the nucleus and the DM mass. For a 10 GeV DM mass and a recoil 
energy Enr of a few KeV, we get for the nuclei we have analysed (Ge, Na, I) that Vm ^ km/s, so the 
expansion is rather accurate (as we have checked explicitly for a Maxwellian halo in appendix A). 
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The unit vectors 'Vg{vm) and 'Vcivm) give the corresponding weight averaged DM wind 
directions in the earth's rest frame. In general they point in different directions and can 
depend on Vm- For the Maxwellian halo they point in the same direction and are equal 
to Vg = \g — ~Vsun- We will treat such special cases in the next section. The positive 
functions g{vm) and G{vm) are given by 

p27r pi 

g{vm)^ J dip J dcos'df{vm,i^,v)cos'd , (16) 

p /•27r rl 

G{vm)= dv d(p dcos'd'f{v,'d',ip)cos'd' , (17) 

J Vm Jo J — 1 

where is the angle between v and V(,(vg). The time dependent halo integral is thus 

given by 

5r]{Vm,t) = Ve(t) ■ [^G{Vm)G{Vm) " Vg(t'm)^^mfi'(^^m)] • (18) 

As already mentioned, the form of Ve(t) from eq. (4) implies that our approximations lead to 
strictly sinusoidal modulations, justifying the ansatz in eq. (9), 5rj{vm-i t) = A^{vm) cos 2n{t — 
to). Using eq. (18), the modulation amplitude can be constrained in the following way: 

Arj{Vm) < Ve [^^^(^m) + G{Vm)] ■ (19) 

Note that the two terms in eq. (18) will contribute to the modulation amplitude proportional 
to sinctg and sincKG, respectively, where ag^a is the angle between Vg^o and the direction 
orthogonal to the plane of the earth orbit, i.e., cosa^^G = ^g,G ■ ©3, with 63 = ei x 62, and 
ei,2 given below eq. (4). In general ag^c will depend on Vm- To derive the bound we have 
assumed the maximal possible effect, corresponding to sina^ = sinac = 1. Since we do not 
know the relative sign of the contributions from g{vm) and G{vm) we have to take the sum 
of the moduli. 

The function g{vm) is bounded from above by 

dip dcosi)f{vm,^,'&)\ cost!}] < / d(pdcos'df{vm,<p,'&)^ 3 — , (20) 

J Vm dVm 

where the last equality follows from eq. (12). Note that r]{vm) is a monotonously decreasing 
function and therefore, dfj{vm)/dvm < 0. Using again the last equality above, also a bound 
for G{vm) can be derived: 



Givm)<-I dv~ = — - I dv-r- (21) 



The equality in (21) follows from integration by parts. The inequalities (20) and (21) are 
saturated if /(v) oc i.e., the hypothetical situation that all DM particles have the same 
direction and their velocities have no transversal component. Using eqs. (20) and (21) the 
bound on the modulation amplitude (19) becomes 



dy ^ Vivm) _ f ^^r]{v) 

dVm Vm Jvrr, '"'^ 



(Assumption 1) , (22) 
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where the first two terms on the r.h.s. are positive and the third is negative. If the 
DM scattering rate R is measured, the r.h.s. is fully determined experimentally by 
vi'^m) = R{Enr)/CF'^{Enr), and cau be compared to the observed modulation through 
= Aji/CF^{Enr), see section 5 for details. Note that the phase of the modulation (which 
may vary with Vm) does not appear in the bound (22). The bound applies on the modulation 
amplitude, irrespective of the phase. Eq. (22) is one of the main results of this paper. The 
bound is rather general and holds under the very mild Assumption 1 specified above. For it 
to be saturated the DM velocity distribution at the position of the solar system would need 
to be highly nontrivial. For instance, the bound in (19) can be saturated if there is a DM 
stream aligned with the ecliptic and it is strong enough so that it dominates the velocity 
distribution ai v — Vm- Note that even in this case, for different the bound will not be 
saturated. To saturate in addition eq. (20) or (21), the stream should have no transversal 
velocity dispersion. Any modulation signal which violates, or even just saturates, this bound 
is very unlikely to have a DM origin. 

4. BOUNDS ON THE ANNUAL MODULATION FOR SYMMETRIC HALOS 

Assumption 2: In addition to Assumption 1 we now assume that defined in eq. (14) 
is independent of Vm, i.e., we assume that there is a constant direction Vhaio = Vg governing 
the shape of the DM velocity distribution in the sun's rest frame. 

Prom this assumption it follows immediately that Vq is also constant and equal to Vhaio, 
so that 

9{vm)^-^, (23) 

dVm 

and eq. (18) becomes 

SviVm, t) = -Ve(t) • Vhaio [Vm9iVm) " G{Vm)] ■ (24) 

The crucial difference to the general case (18) is that we were able to pull the velocity vector 
in front of the bracket. The functions Vmg{vm) and G{vm) are both positive. Their relative 
sizes determine whether the bracket is positive or negative. For small the function G{vjn) 
typically dominates and we get an extra half a year shift in the peak of the modulation (for 
the Maxwellian halo the peak would then be in December). For larger Vm the boundary 
term (fm) dominates and the whole bracket is positive (and thus for the Maxwellian 
halo the peak in this case is in June, see appendix A). 

Assumption 2 is fulfilled if /(v) obeys certain symmetry requirements that we can deduce 
from eq. (14). For a given Vm we chose a coordinate system such that = (1, 0, 0), and 

j d\f{v)vy5{v -vJ^O, j d\f{v)vj{v -vJ^O. (25) 

Assumption 2 implies that these relations hold for any Vm in the same coordinate system. 
The distribution in Vx (i.e., in the direction of Vhaio) can be arbitrary, and can for instance 
even include several peaks, as long as Vg does not flip sign. Therefore, we have to require 
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that the integral over the half-sphere with > is larger than the one with < for all 

Vm- 

I d'vf{^^)vj{v -Vm)< I d?vf{^^)\vMv - Vn.) • (26) 
J Vx<0 J v^>0 

One possibility to satisfy the condition (25) is a symmetric velocity distribution, with 
f{v^, Vy, v^) = f{v^, -Vy, t;^) and f{v^, Vy, v^) = f{v^, Vy, -v^) for all v^. Eq. (25) can also be 
satisfied for distributions asymmetric in Vy and/or v^, however, in this case the asymmetry 
has to be such that the cancellation between Vy^z > and < happens for all radii Vm- 

Assumption 2 is fulfilled for the standard Maxwellian halo, as well as for any other 
isotropic velocity distribution. Up to small corrections due to the peculiar velocity of the sun 
it holds also for tri-axial halos, and covers also streams parallel to the motion of the sun, such 
as a dark disc co-rotating with the galactic stellar disc [19]. Note that for all those examples 
the DM direction Vhaio is aligned with the motion of the sun (up to the peculiar velocity 
that leads to sub- leading corrections). Let us introduce this as an additional assumption: 

Assumption 2a: In addition to Assumption 2 we require that the preferred direction 
Vhaio is aligned with the motion of the sun relative to the galaxy. 

As mentioned above, for many realistic cases fulfilling Assumption 2 also this additional 
requirement is fulfilled. An exception would be the situation when the DM density at 
the sun's location is dominated by a single stream from an arbitrary direction and the 
contribution of the static halo is negligible. 

Let us now use eq. (24) to derive a bound on the modulation. We have Ve(t) ■ Vhaio = 
— fe sin ahaio cos(t — to), where to is now independent of Vm- Here tthaio is the angle between 
Vhaio and ei x 62, i.e., the projection of Vhaio on the ecliptic plane, with sinohaio > 0. 
Assumptions 2 and 2a thus imply that the phase of the modulation is independent of Vm 
(and therefore independent of Enr). As discussed above this is up to a sign flip of the square 
bracket in eq. (24) that can happen due to the two competing terms. To take this into 
account we now define 

5r]{Vra, t) = A'^iVrr,) COs(t - to) , (27) 

where to is constant and we allow a sign change for A'^{vm)- This is different from A^, 
which has been defined to be positive in eq. (9). While for Assumption 2 the phase is 
arbitrary but constant. Assumption 2a also gives a prediction for the phase - that the 
maximum (or minimum) of the event rate is around June 2nd. This can be checked in the 
experiment by looking at the annual modulation phase at different energy bins. Hence, from 
the experimental information on the phase we can conclude whether Assumptions 2 or 2a 
may be justified and whether it makes sense to apply the corresponding test. 

A useful bound on the modulation can be obtained by first integrating eq. (24) over Vm, 

/ dvmA' (vm) =Ve sin ahalo dv ) - G{Vm)\ , . 

Jvml Jvml (,2oj 

=feSinahalo [VmlG{Vml) - Vm2G{Vm2)\ , 

where the second equality is obtained by integrating Vmgiym) by parts using givm) = 
—dG/dVm, eq. (23). Note that the A' and ahaio are defined in such a way that A' is 
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positive for large Vmi above the last sign flip. Experimentally, this is at present the most 
interesting region of the parameter space. Both putative modulation signals, at CoGeNT 
and DAMA, have a peak close to June, and no half year phase change is seen within the 
observed energy range. ^ 

In general we do not know which of the two VmG{vm) terms in (28) dominates. Dropping 
the smaller of the two and using the bound (21) we arrive at 



Vm2 



dVmA'^iVm) 



< sincthaio max 



1l(.Vml) - Vml 



(29) 



Assumption 2 allows for arbitrary directions of the DM wind, therefore we need to use 
Isincvhaiol 1 above, while for Assumption 2a one has sincthaio — 0.5. Note that a sign 
flip of A'^ will lead to cancellations in the integral on the l.h.s. of eq. (29) and make the 
bounds weaker. An observation of such a sign flip in the modulation amplitude would be a 
strong experimental evidence that the modulation is due to a DM signal. If such a sign flip 
is observed, one might be able to obtain stronger constraints by applying the bound for the 
region below and above the sign flip separately, to avoid the cancellations. 

At present there is no indication of such a sign flip in which case the bound (29) simplifies 

to 



/ dVmA!^{Vm) < Ve 



V{Vml) - Vml / dv 



ri{v) 



(Assumption 2) , 



(30) 



and 



dVmA'^iVm) < 0.5' 



ViVml) - Vml / dv 



n{v) 



Vml 



(Assumption 2a) . (31) 



The bounds (22), (30), (31) are the central results of this paper. In the following we will 
refer to them as 



Eq. (22) "general bound" 

Eq. (30) "symmetric halo" 

Eq. (31) "symmetric halo, sin a = 0.5" 



(Assumption 1), 
(Assumption 2), 
(Assumption 2a). 



The term "symmetric halo" should be understood in the sense of eqs. (25) and (26). 

For completeness let us also mention an unintegrated bound, even though we will not use 
it for the numerical analysis. Using eq. (24) we have for the modulation amplitude 



Ar,{Vrn) = Ve\ siuahalol \Vm9{Vm) " G{Vr, 



(32) 



Note that the minus sign between the two terms is conserved, while in the general case, 
eq. (19), one is forced to sum the two terms in the bound. From eqs. (20) and (21) one then 
obtains the following bound on the modulation 



Ar^{Vm) < Ve\ siu Qihalol Hiax 



drj 



dvr. 



dv 



r]{v) 



(33) 



^ Note that the lowest energy bin in DAMA shows a somewhat smaller modulation amplitude, which might 
be an indication of a phase shift below the threshold. 
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A similar bound has been obtained recently in [18]. If eq. (B4) of [18] is expanded in u our 
bound (33) is obtained, if the derivative term dominates. However, eq. (B4) of [18] seems 
to neglect the possibility that the second term in eq. (33) dominates. This assumption is 
justified if data show no phase flip in the modulation and correspond to Vm above the last 
phase flip {A'^ > 0). In this case eq. (33) becomes 

A'^ivm) < -Vesina^^io-, — • (34) 

CLVm 

Note that this bound is complementary to the ones from eqs. (30), (31). By taking the inte- 
gral over the modulation amplitude those bounds probe global properties of the amplitude 
over the considered energy range, whereas eq. (34) bounds the local size of the modulation 
amplitude at a given value of Vm- Both bounds are necessary conditions which a modulation 
signal with a DM origin has to fulfill. 



5. APPLYING THE BOUND TO DATA 

In this section we show how the bounds (22), (30), (31) obtained in Vm space for halo 
integrals can be applied to observable quantities, the unmodulated rate R and the ampli- 
tude of the modulation of the rate, Aj^. The halo integral and the DM scattering rate are 
proportional to each other, see eq. (6). The relation is comphcated by the fact that exper- 
iments typically do not observe the recoil energy Enr directly. The nuclear recoil energy 
is related to the observed energy through a quenching factor Q. In the case of CoGeNT 
and DAMA the observed energy E^e is measured in electron equivalent and is related to the 
recoil energy hj Q = Eee/Enr- In general this is a nonlinear equation, as Q also depends 
on the recoil energy. Finally, in an experiment data is reported in bins, and the continuous 
bounds derived above have to be integrated over the bin sizes. We relegate the details of 
the derivation to appendix B, and quote below only the final results. 



5.1. Single target detector 



Let us first assume that the target consists of a single material, as is for instance the case 
for CoGeNT. We denote the modulation amplitude and the unmodulated rate in bin i as Ai 
and Ri, respectively, both in units of counts/day/kg/keVee. The general bound (22) then 
becomes 



Ai < V, 



N 



Ri{ai + Pi) - Ri+ia'i^^ - {k). ^ Rjjj 



3=1 



(Assumption 1) 



(35) 



and the bounds for the symmetric halo are 



TV 



Aj Xj < f e sin a 



j=i 



N 



RiVi - {Vm)i ^ Rjlj 



(Assumptions 2, 2a) 



(36) 



where the bounds (30) (Assumption 2) and (31) (Assumption 2a) are obtained for sin a = 1 
and 0.5, respectively. The bin index i runs from 1 to A^, while the rates Ri are zero for 
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i > N. The coefficients a^, a', f3i, (k)- , 7^, a;,, yi, {vm)i are given in eqs. (B9) and (B12). They 
are known quantities calculable in terms of the form factor F{Enr), quenching factor Q{Enr), 
and Jacobians needed for changing the variable from to Eg^. They depend on the DM 
mass m^^ via the reduced mass /i^-a, eq. (2), needed to convert into Enr- The dependence 
on the scattering cross section and the local DM density is encoded in the constant factor 
C, eq. (6). This factor cancels completely, as expected, since it is a common factor for the 
modulation as well as for the rate. 

The rates Ri in the bounds (35) and (36) are the unmodulated scattering rate induced 
by DM without including any backgrounds. In a "background free" experiment, where the 
full observed event rate is due to DM, the bounds can be apphed as they are. Here we want 
to be more conservative, and consider also the situation where an unknown background 
may contribute to the unmodulated rate. The remaining assumption is then just that the 
background itself is not modulated, only the DM signal. In each bin a fraction cui of the 
observed count rate TZi is due to DM, i.e.. 



Ri = TZiUJi , <uji <1, 
so that we can replace Ri — >■ TZiUJi in eqs. (35) and (36) to obtain. 



(37) 



Ai < Ve. 



N 



j=i 



N 



AjXj < Ve sin a 



j=i 



N 



( Assumption 1) , (38) 
(Assumptions 2, 2a) . (39) 



Now we have to find a set of cjj (i = 1, . . ., A^), such that the bound becomes the "weakest". 
In the following we describe two different procedures for this task. Procedure 1 is easy to 
implement but gives slightly weaker bounds, whereas procedure 2 involves an optimization 
algorithm but gives somewhat more stringent bounds. 

Procedure 1: We can sum eq. (38) from bin i to N and drop the last term in the square 
bracket, since {k.)^ and 7^ are positive: 



N 



N 



N 



(Assumption 1) , (40) 



where the last inequality was obtained by setting a;^ = 1. Similarly we can drop the second 
term from eq. (39) and obtain 



N 



y^ AjXj < Ve sin a TZiUiyi < Ve sin a TZiyi (Assumptions 2, 2a) , 



(41) 



j=i 



where the last inequality holds for yi > These bounds have to be satisfied for all bins i. 



According to eq. (B12) yi can also become negative. In that case the bound (41) is always violated for 
any coi > 0. 
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Procedure 2: A somewhat stronger bound can be obtained by searching for the optimal 
choice for the Ui, without dropping the last terms in eqs. (38) and (39). We present here 
a method based on a least-square minimization (an alternative procedure is outlined in 
appendix C). Let us denote the r.h.s. of eqs. (38) and (39) as Bi which are functions of cUj. 
Then we can construct from eq. (38) the following least-square function 

^ //I _ R \ 2 

X^ = Y,\ ' A ' ) ^{^^- B^) (Assumption 1) , (42) 

i=l V I / 

where af" is the la error on Ai (errors on Bi are typically much smaller and we neglect them 
here). This can now be minimized with respect to Uj under the condition < < 1. 
The G function takes into account that there is only a contribution to if the bound is 
violated. Hence, X"^ will be zero if the bound is satisfied for all bins. A non-zero value of 
indicates that the bound is violated for some bin(s), weighted by the corresponding error 
in the usual way. In the case of eq. (39) one has to take into account that the modulation 
amplitude in each bin is used several times, leading to correlated errors for the l.h.s.: 

N 

= ^ {Ai - Bi)Sr.\Aj - Bj)Q{Ai - Bi)Q{Aj - Bj) (Assumption 2,2a) , (43) 



where 



N N , , , , N 



k=i k=l fc=max(i,j) 

While non-zero values of the functions (42) and (43) can be considered as a qualitative 
measure for the violation of the bound, the precise distribution of them should be determined 
by Monte Carlo studies. From the definition one can expect, however, that they will be 
approximately distributed if the bound is violated. 



5.1.1. CoGeNT 

Let us consider now the modulation signal reported by the CoGeNT experiment at 
2.8(7 [4]. It has been pointed out that for specific assumptions on the halo (standard 
Maxwellian) there is a tension between the modulated and unmodulated rate in Co- 
GeNT [20, 22, 23]. Recent analyses on the CoGeNT modulation signal can be found in 
refs. [24, 25], see also [18]. Here we use the CoGeNT data for a case study and apply the 
above bounds assuming spin-independent elastic scattering. For the germanium quenching 
factor we use £'ee[keV] = 0.199(i?„r[keV])^'^^ [26]. We adopt the Helm parameterization 
of the spin-independent form factor, F{Enr) = 3e~'' * /2[sin(gr) — gr cos(gr)]/(gr)^, with 

= 2mAEnr and s = 1 fm, r = y/R^ — Ss^, R = 1.2^1^/'^ fm. We use the data from fig. 6 of 
ref. [20] where the total rate TZ and the modulation amplitude are given in four bins of 
Egg between 0.5 and 3.1 keV. We reproduce the data in table I.^ 



^ We thank Mariangela Lisanti for providing us the data from fig. 6 of ref. [20] . 
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F, bins fkpVl 


MnH ("Ass 1 2) 




TTTTmnH ra'I'P 




0.5-0.9 


1.41 ±0.79 


0.90 ±0.72 


12.33 ±0.52 


5.29 ±0.52 


0.9-1.5 


0.84 ±0.59 


0.37 ±0.55 


4.33 ±0.39 


3.36 ±0.39 


1.5-2.3 


0.46 ±0.24 


0.48 ±0.22 


2.76 ±0.16 


2.76 ±0.16 


2.3-3.1 


0.66 ±0.24 


0.27 ±0.23 


2.83 ±0.17 


2.83 ±0.17 



TABLE I: CoGeNT data on modulation amplitude and unmodulated rate [4] in cnts/day/kg/keV, as 
reported in . 6 of ref. [20]. The modulation has been extracted allowing for indpendent phases in each bin 
(Ass. 1) and for a constant but arbitrary phase (Ass. 2) (which lead to very similar amplitudes), and by 
requiring the maximum at June 2nd (Ass. 2a). In the last column we show the preliminary surface events 
corrected unmodulated rate [21]. 




8.5-Q.9 0.9-1.5 1.5-2.3 2.3-3.1 

fee (KeV) 



FIG. 1: Procedure 1 upper bound compared to the integrated modulation amplitude from the CoGeNT 
data, from fig. 6 of [20]. The red dots correspond to the l.h.s. of eq. (40). The black diamonds correspond 
to the upper bound obtained under Assumption 1, r.h.s of eq. (40), for DM masses of rrij^ = 7, 10, 20 GeV 
(from bottom to top) . The bins on the horizontal axis indicate the bin i from which we start to sum the 
data. Error bars correspond to la. 

The modulation amplitude has been extracted in three different ways in [20]. First, 
by fitting independently the modulation phase for each energy bin, second, by assuming a 
constant phase for all bins, and third, by fixing the phase such that the modulation maximum 
is on June 2nd. These are precisely the requirements corresponding to our Assumptions 1, 
2, 2a, respectively. We can thus use the appropriate data on the modulation for each of the 
three assumptions. The modulation amplitudes in the first and the second case are very 
similar, while in the third case (forcing the phase to equal June 2nd) the amplitudes are 
lower and error bars are larger. 

Choosing three values of the DM mass as examples, we show in fig. 1 the bound on the 
integrated modulation amplitude for Assumption 1, and using procedure 1. The bin labels 
on the horizontal axis give the iih energy bin, which is the lower limit of summation in 
eq. (40). The data on amplitudes are shown as red dots, while the bounds are shown in 
black. Whenever a red dot lies above one of the bounds (within errors), the DM hypothesis 
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FIG. 2: Bounds on the CoGeNT modulation amplitude for the three assumptions about the DM halo (1, 
2, 2a), and using eqs. (40) and (41) (procedure 1). We show the l.h.s. (integrated modulation amplitude) 
minus the r.h.s. (upper bound) divided by the error on the l.h.s., as a function of DM mass. We sum the 
data starting from bin i = 3,2 and 3 for Assumptions 1, 2 and 2a, respectively. The bounds are violated in 
the regions where the curves are above zero, which are shaded in the plot. 



is disfavoured. In fig. 1 this happens for the case of DM mass = 7 GeV. Such light DM 
thus cannot be the source of the modulation, even under the very general assumption on 
the DM halo adopted here. 

We see from fig. 1 that the strongest restriction comes from the bins z = 2 or 3. In 
fig. 2 we show the constraints on DM mass that follow from the bounds on the modulation 
amplitudes in these two bins. The bounds are obtained from procedure 1 for the three 
Assumptions 1, 2, 2a. We observe that the bound becomes stronger for smaller DM masses. 
This behaviour can be understood from how the coefficients defined in appendix B depend 
on the reduced mass fj,^A- The bounds on amplitudes are violated for DM masses below 
10, 27, 33 GeV, respectively, for the three assumptions. The lower bounds on are 
summarized in table II (left part), where we also give the bounds at la and 2a. Notice that 
due to the smaller modulation amplitudes and larger errors when extracted from the data 
under Assumption 2a, at 95% CL the corresponding bound becomes weaker than the one 
from Assumption 2 and equal to the one from Assumption 1, although naively one would 
expect the opposite. This can be traced back to the fact that, under Assumption 2a, the 
modulation phase is forced to take the value of June 2nd which is not the one preferred by 
the data. The extracted modulation signal then gets weaker and consequently the bounds 
are more easily satisfied. We have also checked that bounds for Assumptions 2, 2a derived 
from eq. (34) give always weaker limits than the ones discussed here, which are based on 
eq. (30) and eq. (31). 

In fig. 3 we show the functions according to procedure 2 as a function of the DM mass. 
For a given value of we minimize X"^ numerically with respect to the Uj. The conclusion 
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unmodulated rate from [4] 


corrected unmod. rate [21] 




Mean 


68% 


95% 


X2 < 1 


Mean 


68% 


95% 


X2 < 1 


Ass. 1: general bound 




8.5 


6 


3 


7.3 


10 


6.5 


3 


10 


Ass. 2: symmetric halo 




24 


14 


6 


18 


43 


25 


12.5 


37 


Ass. 2a: sym. halo, sin a 


= 0.5 


27.5 


13.5 


3.5 


16 


59.5 


23 


3 


35 



TABLE II: Lower bounds on the DM mass in GeV, from the requirement that the modulation amplitude 
in CoGeNT is consistent with the upper bound from the unmodulated rate, according to the Assumptions 
1, 2, 2a on the DM distribution. The bounds are obtained from procedure 1, requiring that eqs. (40) or (41) 
are satisfied for the mean value, or the 68% and 95% CL limits. The bounds labeled "X^ < 1" are obtained 
from procedure 2 by requiring that X'^ defined in eqs. (42) or (43) is less than 1. In the left part of the table 
we use the published unmodulated event rates from [4] , whereas for the right part of the table we adopt the 
preliminary results on surface events contamination at low energies from [21]. 
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FIG. 3: Bounds on the CoGeNT modulation amplitude for the three assumptions about the DM halo (1, 
2, 2a), always using procedure 2. We show defined in eqs. (42) and (43) as a function of DM mass. The 
thin dashed curve is for illustrative purpose only; it corresponds to the data and errors on the modulation 
amplitude extracted with a free phase (Assumption 2) but using the bound for Assumption 2a, see text for 
details. If we assume that X^ is distributed as a with 1 d.o.f. (see below) X^ — 1 (4) corresponds to 68% 
(95%) CL. 

is similar to procedure 1, leading to similar lower bounds on the DM mass. Requiring that 
< 1 one finds > 7.3, 18, 16 GeV for Assumptions 1, 2, 2a, respectively. We observe 
again the unusual situation that Assumption 2a leads to a weaker constraint, due to the 
less significant signal for the modulation. To illustrate this effect we show in fig. 3 with 
a thin-dashed curve also X'^ that would follow from the hypothetical situation where the 
modulation amplitude would be as strong as in the case of Assumption 2. That is, we allow 
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FIG. 4: Bounds on the CoGcNT modulation amplitude for the preliminary surface events corrected 
unmodulated rate [21] for the three assumptions about the DM halo (1, 2, 2a). The left plot shows the 
bounds according to procedure 1, whereas the right plot shows defined in eqs. (42) and (43) (procedure 2). 
The thin curves in the right panel are the bounds without surface event subtruction reproduced from fig. 3 
for the purpose of comparison with the thick curves obtain with the surface event corrected rates. 

for an arbitrary energy independent phase, but assume that this phase turned out to be 
June 2nd (in contrast to the real CoGeNT data), so that we could apply the bound from 
Assumption 2a to these data. We see that in this case one would obtain a much stronger 
bound, disfavoring DM masses up to 60 GeV. This example shows the potential of our 
method in cases of a strong modulation signal in the data. If the signal for the modulation 
itself is weak, the bounds we derived will also give only weak constraints. 

A recent re-analysis of CoGeNT data indicates that a significant fraction of the event 
excess at low energies could be due to surface events [21]. Assuming that surface events are 
not modulated, the unmodulated rate will be reduced after subtructing the surface events, 
while the modulation signal will remain, leading to a strengthening of our bounds. Here we 
estimate this effect by using the preliminary result for a surface events rejection efficiency 
shown on slide 19 of the presentation in [21] (red curve). Averaging this curve for the bins 
used in our analysis we find that the unmodulated rate in the first and second bins are 
reduced by a factor 0.43 and 0.78, respectively, while the other two bins are not effected. 
This leads to the reduced event rates shown in the last column of tab. I. Fig. 4 shows 
the plots equivalent to figs. 2 and 3 but using the surface events corrected unmodulated 
rate. The general bound remains essentially the same. In this case the strongest limit (from 
procedure 1) for the uncorrected rate comes from summing data starting from bin 3, whereas 
in the surface events corrected case it comes from bin 2, with only a minor reduction (79%) of 
the rate, leading to a very similar limit. In contrast, for Assumptions 2 and 2a significantly 
stronger bounds are obtained, with the limit coming now from summing from bin 1, with a 
43% reduced rate. The surface events corrected bounds are summarized in the right part of 
table. II. 

Due to the non-standard definition of the X"^ functions (42) and (43), involving the G- 
function, the actual distribution of them is not clear a priori. Therefore we have performed 
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FIG. 5: The probability to obtain a larger than the one obtained from CoGeNT data, as determined 
by Monte Carlo simulation. Dashed curves correspond to surface event corrected data. 



a Monte Carlo study in order to determine the distribution. For a given DM mass vi^e 
first determine tlie optimal set of by minimizing the X'^ . For those Ui vi^e assume that 
the bound is saturated and we simulate a large number of pseudo-data taking the r.h.s. of 
eqs. (38), (39) as mean value for a Gaussian with standard deviation given by the actual 
error on the modulation amplitude.^ For each random data set we calculate the value 
and obtain therefore the distribution of assuming that the bound is saturated. Then we 
can compare the of the actual data and calculate the probability of obtaining a larger 
than the observed one. Note that this procedure is conservative, since we assume that the 
mean value of the random data is the bound itself. If the "true" mean values are smaller 
than the bound, the distribution of X^ would be shifted towards smaller values, leading to 
more stringent bounds. 

We show the results of this calculation in fig. 5. The probabilities obtained in this way 
are in qualitative agreement with the numbers reported in tab. II. Under assumptions 2 and 
2a, surface event corrected data is inconsistent with the DM hypothesis for any DM mass 
at the 97% and 90% CL, respectively. We find that if the bound is violated (i.e., for small 
DM masses) the X^ distribution is close to a x^-distribution with 1 d.o.f.. Large deviations 
from a x^-distribution occur if X^ is close to zero (i.e., for larger DM masses). In this case 
the X^ distribution is strongly peaked at small values close to zero. This indicates that in 
such situations it is in most cases possible to find a set of a;, such that X^ becomes zero. 
This is the reason why non-trivial constraints are obtained under assumptions 2, 2a for 



^ In the case of Assumptions 2 and 2a we proceed iteratively. Starting from eq. (39) for i = N we generate 
Aj\j, and then simulate successively Ai^i- This is necessary to obtain the correct random properties of 
the l.h.s. of eq. (39). 
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surface events corrected data at large masses, although the values of X"^ are relatively small, 
compare fig. 4. 



5.2. Bounds for multi-target experiments 

So far we have restricted the discussion to the situation where only one target nucleus is 
present, as for example Ge in CoGeNT. Let us now generalize our bounds to experiments 
where several elements are used as DM target. In this case both the modulation amplitude 
as well as the unmodulated rate will receive contributions from each clement: Ai — 
and Ri = i?", where n labels the different target elements and i energy bins. The 
bounds from procedure 1 hold for each of the elements separately, where in general all the 
coefficients appearing in the equations will depend on the nucleus type. Summing eqs. (40) 
and (41) over n and using the fact that for positive Oj and bi, the following inequality holds, 
Zli^i^i < {J2i(^i)(^jbj), we derive: 



N 



3=1 

N 



N 



(Assumption 1) , 



Aj min„ (x^) < Ve sin a TZi ||/"| (Assumptions 2, 2a) . 



(45) 



(46) 



On the l.h.s. of the last inequality it is necessary to take the minimum of between all 
the nuclei present in the target, for each bin. We used that and Xj are positive, see 
appendix B. 

It is possible to generalize also the X"^ method to the case of multiple targets. For each 
nucleus the bounds (38) and (39) apply separately We define the amplitude from element 
n in the bin i as A^ = e'^Ai, and similarly for the rate i?" = u^TZi. Then one can construct 
a X^ similar to the one in the previous section, from l.h.s. minus r.h.s. of the following 
inequality: 



e'lAi < Ve 



N 



(Assumption 1) . (47) 



This X"^ is minimized over a;" > 0, > given the constraints J2n '^i — ^ ^i=^ 
for each i. A similar relation can be set up for Assumption 2. 



5.2.1. DAM A 

Let us now apply above bounds to the DAMA result, based on data from a Nal detector. 
DAM A observes a relatively large unmodulated count rate of about 1 cnt/kg/day/keV com- 
pared to the modulation amplitude of about 0.02 cnt/kg/day/kcV. In contrast, CoGeNT 
reports a modulation amplitude in the range of 0.5 to 1 cnt/kg/day/keV, compared to an 
unmodulated rate of 3 to 4 cnt/kg/day/keV. Because of this much smaller modulation am- 
phtude compared to the unmodulated rate for DAMA we expect the signal to be consistent 
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FIG. 6: Upper bound compared to the integrated modulation amplitude for DAM A data for Assumption 1 
(left) and Assumptions 2 and 2a (right). We assume quenching factors q = 0.3 and q = 0.09 for Na and 
I, respectively. The red dots correspond to the l.h.s. of eqs. (45) and (46), and the purple triangles to the 
r.h.s.. In the left panel we show the bounds for DM masses of = 3,20,50 GeV (from bottom), in the 
right panel for = 3 GeV. Bins are 0.5 keV wide and we sum all bins starting from bin i shown on the 
horizontal axis up to 7 keVee. Error bars are negligible and are not shown for clarity. All dark matter 
masses are compatible with the modulation. 

with our bounds. It is also well known that a consistent fit of the DAMA modulated and un- 
modulated data is possible e.g., for a Maxwellian halo, although in certain parameter regions 
the unmodulated rate does provide a non-trivial constraint, see e.g. [8, 11]. In fig. 6 we show 
the bounds for DAMA, assuming quenching factors of 0.3 for Na and 0.09 for I (constant in 
energy) using the spectral data on the modulated and unmodulated rates from [3] . We find 
that both the general bound, eq. (45), as well as the bounds for a symmetric halo, eq. (46), 
are consistent with the modulation even for a DM mass as small as 3 GeV. As noted above 
the bounds get weaker for larger DM masses, hence they do not place any relevant constraint 
on the dark matter interpretation of the DAMA modulation signal. 

In deriving the multi-target bounds eqs. (45) and (46) we had to use inequalities related to 
the presence of several contributions to the rates, which makes the bounds somewhat weaker 
than single target bounds. In case of Nal we have a relatively large hierarchy between the 
two targets since A = 23 for Na and A = 127 for I. Therefore, one may consider also 
the situation that only one of them dominates the signal. If scattering on both targets 
is kinematically allowed, iodine will dominate because of the enhancement of the spin- 
independent scattering cross section (which we will assume here). For low DM masses, 
scattering on iodine may be kinematically forbidden since is larger than the escape 
velocity of the halo and in that case the signal is provided only by scattering off Na (this 
requires additional assumptions on the escape velocity, which in general are not necessary 
for our bounds). We show the bounds for Assumptions 2 and 2a assuming that scattering on 
either Na or I dominates in fig. 7 for some representative DM masses. In both cases only very 
weak constraints are obtained. While scattering on Na is consistent down to ^ 3 GeV 
for Assumption 2a, scattering on iodine is consistent for > 30 GeV under Assumption 2. 
(Note that below around 30 GeV typically scattering on Na dominates.) 

As expected, we conclude that based solely on our bounds the DAMA signal is compatible 
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FIG. 7: Upper bound compared to the integrated modulation amplitude from DAMA data for Assump- 
tions 2 and 2a, assuming that DM scatters only on sodium (left) or iodine (right). The red dots correspond 
to the l.h.s. of eq. (41) and the purple triangles to the r.h.s., for a DM mass of — 7 GeV (left) and 
30 GeV (right). Bins are 0.5 keV wide and we sum all bins starting from bin i shown on the horizontal axis 
up to 7 keVee. Error bars are negligible and are not shown for clarity. 

with assuming dark matter scattering being its origin. We do not expect that this conclusion 
will change significantly if bounds from procedure 2 according to eq. (47) were applied. 
Therefore, we limited ourselves to the discussion based on procedure 1 for DAMA. 



6. SUMMARY AND DISCUSSION 

The annual modulation is probably the most distinctive signature of dark matter and is 
playing (and will certainly play) a central role in revealing its existence and nature. However, 
it is not always clear that the modulation detected by a DM direct detection experiment 
is caused by the dark matter particles and not by any other unknown background source. 
In this work we have presented a consistency check for the amplitude of the modulation 
and the unmodulated count rate, which is a necessary (but not sufficient) condition for DM 
being the origin of an observed annual modulation. 

We have derived upper bounds on the energy integrated modulation amplitude in terms 
of the unmodulated rate by expanding the DM velocity integral to first order in the earth's 
velocity Ve ~ 30 km/s. It holds for a wide class of particle physics models where the 
differential scattering cross section da/dEnr is proportional to l/v"^. Although we have only 
focused on elastic scattering the method can also be generalized to the inelastic case. The 
important aspect of our work is that our bounds hold for very general assumptions about 
the DM velocity distribution /(v) in the sun's vicinity. We have presented bounds under 
the hypothesis of a single DM species and the following assumptions on the DM halo: 

• Assumption 1, bound in eq. (22): We assume that the only time dependence is 
induced by the rotation of the earth around the sun. The halo itself is static on the 
time scale of months to years and spatially constant at the scale of the sun-earth 
distance. Otherwise the DM velocity distribution can have an arbitrary structure, 
including, for example, several streams coming from various directions. In order to 
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saturate this bound the halo has to be very pecuhar, with rather unreahstic properties. 

• Assumption 2, bound in eq. (30): In addition to Assumption 1 we assume that 
there is just one preferred direction of the DM velocity distribution in the rest frame 
of the sun (independent of the minimal velocity Vm in the halo integral). This requires 
certain symmetries of the velocity distribution, which are specified in eqs. (25) and (26). 
It covers typically situations where the DM distribution is dominated by one single 
component, which may come from an arbitrary (but constant) direction. Assumption 2 
requires that the phase of the modulation is independent of the recoil energy (up to a 
phase flip of half a year) . 

• Assumption 2a, bound in eq. (31): We require that the preferred direction from 
Assumption 2 is ahgned with the motion of the sun. This is fulfilled for any isotropic 
halo, and also for tri-axial halos up to corrections due to the peculiar velocity of the 
sun. Furthermore, it includes the possibility of streams parallel to the motion of the 
sun, such as a possible dark disk. Assumption 2a requires that the maximum (or 
minimum) of the event rate is around June 2nd. 

The theoretical bounds are obtained in terms of the minimal velocity Vm and have to 
be related to observable quantities by translating into recoil energy. We have outlined two 
possible procedures for this task in section 5 for the case of elastic scattering, taking into 
account statistical errors and the possibility that an unknown background may contribute to 
the unmodulated rate, but that it has no time-dependence. As an example we have applied 
the proposed consistency checks to the annual modulation signals reported by the CoGeNT 
and DAMA experiments. While DAMA data are compatible with a dark matter origin for its 
modulation, severe restrictions on the dark matter mass can be set for the case of CoGeNT. 
Applying our bounds we find that the CoGeNT modulation amplitude can be consistent 
with the unmodulated rate at the 68% CL only for DM masses > 7.3, 18, 16 GeV 
for the Assumptions 1, 2, 2 a, respectively. If preliminary results on a possible surface 
events contamination of the unmodulated rate at low energies in CoGeNT are confirmed, 
those bounds would become even more restrictive. In this case, CoGeNT modulation data 
would be inconsistent with the DM hypothesis under assumptions 2 and 2a at about 97% 
and 90% CL, respectively. DAMA has a relatively large unmodulated count rate of about 
1 cnt/kg/day/keV compared to the modulation amplitude of about 0.02 cnt/kg/day/keV. 
Therefore, our bounds are not very stringent and the modulation amplitude is consistent with 
the unmodulated rate. In contrast, CoGeNT reports a modulation amplitude in the range 
of 0.5 to 1 cnt/kg/day/keV, compared to an unmodulated rate of 3 to 4 cnt/kg/day/keV. 
Because of this relatively large ratio between modulated and unmodulated rate our method 
provides stringent constraints in the case of CoGeNT. Several comments are in order: 

(i) In deriving the bounds we assume a certain smoothness of the halo, since we are 
expanding in Ve- Spikes in the velocity distribution much narrower than 30 km/s are not 
covered by our procedure. 

(ii) Our bounds assume that a modulation signal is present in the data. If the sig- 
nificance of the modulation is weak, the bounds are more easily satisfied. This effect is 
explicitly illustrated in sec. 5 by using CoGeNT data under different assumptions regarding 



21 



the modulation signal. The method discussed here cannot be used to establish the presence 
of a modulation, it can only test whether a given modulation signal is consistent with the 
unmodulated rate. 

[iii) In this work we have used the relation between the minimal velocity Vjn and nuclear 
recoil energy E^r implied by elastic scattering, see eq. (2) and appendix B. The bounds 
can be generalised in a straightforward way also to inelastic scattering. However, in that 
case Vm = {rriAEnr/ fJ'xA + 5) / \/2mAEnr has a minimum in Enr-, and therefore there is no 
one-to-one correspondence between Vm and E^r- When translating from Vm to Enr one has 
to take into account that different disconnected regions in Enr can contribute to a given 
interval in Vm- 

(iv) The type of DM-nucleus interaction (i.e., the particle physics) has to be specified 
before applying our bounds. Apart from the dependence of the differential cross section 
da/dEnr, the E^r dependence of the interaction is encoded in the form factor F{Enr)- It can 
describe conventional spin-independent or spin-dependent interactions, but also a possible 
non-trivial g'^-dependence of the DM-nucleus interaction, all of which can be absorbed into 

F{Enr)- 

To conclude, we have presented a consistency check for the amplitude of a DM induced 
annual modulation compared to the unmodulated event rate in a given DM direct detection 
experiment. Our bounds rely only on very mild and realistic assumptions about the DM halo. 
We believe that the proposed method is a useful check which any annual modulation signal 
should pass. A violation of our bounds suggests a non-DM origin of the annual modulation, 
or requires rather exotic properties of the DM distribution, for example very sharp spikes and 
edges. Such extreme features of the halo should have other observable consequences, such 
as surprising spectral features or strong energy dependence of the modulation phase. Such 
features could be used as additional diagnostics, beyond the signatures explored here, which 
are restricted to the energy integrated modulation amplitude, irrespective of the phase. 
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Appendix A: Maxwellian halo 

As a cross check we apply our formalism to the commonly used Maxwellian velocity 
distribution, with a cut off at an escape velocity v^sc- In the galactic rest frame it is given 

by 

/gai(v) oc [exp(-v7'y^) - ex.p{vl^Jv'^)]e{ve^ - v) , (Al) 
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FIG. 8: Comparison for the Maxwehian halo of the exact (solid line) and expanded (dashed) modulated 
and unmodulated halo integrals A^^ and ry, respectively. The modulation flips the phase by half a year below 
the zero of the amplitude. In the lower panel the general bound (22) is also shown (dashed-dotted curve). 
The upper panel shows the relative accuracy defined as (exact - approx) /exact. 

where we use v = 220kms~^ and fesc = 550 km s^^. This distribution is boosted to the sun's 
and earth's rest frames as described in section 2. For the above velocity distribution the halo 
integral ri{vm,t) is known analytically. We define the "exact" modulated and unmodulated 
halo integrals as 



V (Vm) = 2 Iviv^m, Q + r]{Vm, U + 0.5)] , 



A J exact / 



(A2) 



where for the (exact) analytic expressions ri{vm,t) we are are using the results given in the 



exact / 



I A 



/ exact f 



appendix of [27]. Above, is June 2nd and A. 

We first check the accuracy of the Ve expansion. Expanding to zeroth order in Ve gives the 
approximate unmodulated rate (12), while the modulation amplitude is given by expanding 
to linear order in Ve, (32), with sinahaio = 0.5. The comparison with the "exact" expressions 
is shown in fig. 8. Both for the modulation amplitude and the unmodulated rate we find 
excellent agreement, with the differences hardly visible on logarithmic scale. As seen from 
the upper panel, for large regions of Vm the agreement is within few %, apart from the zero 
of the modulation amplitude. Minor deviations appear for large Vm values, where non-linear 
effects due to the cut-off at the escape velocity become important. We see that our expression 
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FIG. 9: The bounds on the modulation amplitude for a Maxwellian velocity distribution and velocity 
dispersions v = 220 km/s (40 km/s) in the left (right) panels. The shaded region shows the integral of 
^ A'^°^'^'^^. We show the bound for a "symmetric halo" from eq. (30), and the bound from eq. (31) which 
assumes that the DM wind is parallel to the motion of the sun ("sym. halo, sin a = 0.5"). The solid curve 
("partial bound") shows the bound eq. (30), but without the second term in the square bracket. 

for captures accurately the phase flip of the modulation. At low Vm the G{vm) term in 
eq. (32) dominates, leading to the maximum of the count rate in December, the modulation 
amplitude becomes zero when the two terms cancel, and the maximum in June at large Vm 
is obtained from the g{vm) term. 

We compare next the bounds with the actual modulation amplitude. The dash-dotted 
(green) curve in fig. 8 shows the general bound, eq. (22), which follows from Assumption 1. 
Clearly, the Maxwellian halo is far from saturating this bound. Much tighter constraints are 
obtained from Assumptions 2 and 2a. Fig. 9 shows the bounds (30) and (31) on integrals 
of modulation amplitudes. The shaded region shows the l.h.s. of the bounds, f dvA', 
where the upper boundary of the integration is chosen so high that the DM signal vanishes. 
Note that A'^ flips the sign at Vm ~ 200 km/s, which explains the maximum of the integral 
around 200 km/s. The curve labeled "partial bound" corresponds to the bound following 
from Assumption 2, eq. (30), but without the second term in the square bracket. We find 
that it is only slightly more stringent than the general bound that follows from Assumption 
1, (22), especially for large Vm- Although it does not capture the behaviour at low Vm it is 
rather similar to the full bound, eq. (30), (dashed curve) for Vm ^ 400 km/s. The bound 
from Assumption 2a, eq. (31) (dash-dotted curve), which assumes sinahaio = 0.5, is roughly 
a factor 2 larger than the prediction for small Vmi and approaches the true modulation 
around Vmi ^ 400 km/s. In the right panel we show that for a Maxwellian distribution with 
very small dispersion {v = 40 km/s) the bound for Assumption 2a is nearly saturated. As 
mentioned in sec. 3, the inequalities used to derive the bound (31) become equalities if the 
velocity distribution /(v) oc 6{i)), which is approximated by the Maxwellian with small v. 
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Appendix B: Translating bounds in Vm space to observable quantities 



In sections 3 and 4 wc have derived bounds involving the halo integral fj{vm) and the 
modulation amplitude Arj{vm)- Here we give details on how to translate these bounds into 
bounds that involve observable rates as functions of nuclear recoil energy E^r- The nuclear 
energy is measured in electron-equivalents, E^e, which is related to Enr through a quenching 
factor Q. We also define a differential quenching factor q, 



where Q 



E. 



Er, 



which then gives the translation between Vm and E^e, 

1 dVrr,. 



ImAEee 



dvjji dEfij. 



C{Eee) dEee dE^r dE^e '^jJ'xA V 2£'ee Q 



rriAQ 1 



(Bl) 



(B2) 



The quenching factor Q is in general energy dependent. If it is energy independent, then 
q — Q, while if data are reported directly in Enr then furthermore q — Q — 1. Note that 
(l,Q,Vm,Cj are all positive. 

The modulation amplitude and rate in units of counts/day/kg/keVee are given as 

CF\EJQ) 



AR{E,e) 



R{Eee) = Hr]{Vm) with K{Eee) 



q{Eee) 



(B3) 



with the constant C defined in (6). To keep the notation simple we use the same symbols 
for rate and modulation whether they depend on recoil energy in keVge or keVn,.- The units 
are denoted in the argument {E^e versus Enr) and should be clear from the context. 

Let us consider first the general bound following from Assumption 1. Multiplying eq. (22) 
by K and converting the integral over Vm into Ege by using eq. (B2) one finds 



AniEee) < Ve 



djRQ 

dEe.e. 



+ R 



d^ 

dE,, 



TUAEee H dE^e 



K 



dE, 



R ^4aQ 

TUAEee 



Eq. (34) for Assumptions 2, 2a becomes similarly: 



(B4) 



A-R{Eee) < Ve Siu QJhalo 



dEe 



+ R 



dC, C, dn 



dEe 



+ 



K dEe 



(B5) 



Note that the constant C containing the cross section and local DM density drops out, 
as expected, since it is a common pre-factor for rate as well as modulation. In practise Ar 
and R are given in bins in Eee- Let us define the bin average of a quantity X{Eee) as 



AE, 



dEeeX{Eee) , 



(B6) 



where Eii and Ei2 are the boundaries of bin i and AEi = Ei2 — Ei\. The observed modulation 
and rate in bin i are then Ai = (Ar)- and Ri = {R)^, respectively. Now we have to perform 
the bin average of eq. (B4). The first term in the square bracket gives 



R{Ea)aE^l) - R{E,2)i{E,2) Rr {Oi ' Rr+i iO. 



+1 



AE, 



AE, 



(B7) 
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For the other two terms in the square bracket we approximate {XY) fa (X) (Y) and the 
integral becomes a sum over bins, up to the last reported bin i — N. Finally we obtain 



Ai < Ve. 



N 



Ri{ai + A) - i?i+iQ;-+i - (k), Y1 Rffj 



j=i 



(B8) 



with 



For eq. (B5) we get 



dE,, 



rriAEee K dEe 



, 7i = 



/ d^ ^ dn 

Ai < Ve sin fthaio [Ri{ai + /3'^ - Ri+ia[^^] with /3[ = ( + - 



dE, 



mAEeein/ ^ 

(B9) 



(BIO) 



Terms with indices larger than N are dropped. We see from the definitions that Q;j,Q;^,7i, 
and (k). are positive. 

A similar calculation for the bounds (30) and (31), following from the Assumptions 2 and 
2 a, leads to 



N 



AjXj < Ve sin a 



N 



RiVi - {Vm)i 



(Bll) 



with 



X, = AE, 



irriAEe 



(B12) 



and sin a = 1 (0.5) corresponds to Assumption 2 (2a). 

Any binning procedure will lead to inaccuracies. We have estimated binning effects by 
slightly changing the prescription. Eq. (B4) can be equivalently written as 



AniEee) < Ve 



^ dEpp 



'2/i2^Q ^ e dK 



rUAEee K dEe 



dEp 



R '^i^IaQ 

''^K rUAEee 



(B13) 



The first term in the square bracket can be estimated as 



dR 
dEpp 



Ri — Ri 



+1 



AEi 



(B14) 



and the other terms are bin-averaged similar as above. We have calculated the CoGeNT 
bounds also using this binning procedure and obtain similar results as with our default 
method. This confirms that inaccuracies due to binning are acceptable for present CoGeNT 
data. 
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unmod. rate from [4] 


corrected rate [21] 




Mean 


68% 


95% 


Mean 


68% 


95% 


Ass. 1: general bound 




10 


5 




12.5 


5 




Ass. 2: symmetric halo 




29.5 


17 


7 


63 


34.5 


16 


Ass. 2a: sym. halo, sin a 


= 0.5 


37.5 


14 




94.5 


30.5 





TABLE III: Lower bounds on the DM mass in GeV, from the requirement that the CoGeNT modulation 
amplitude is consistent with the upper bound from the unmodulated rate, according to the Assumptions 1, 
2, 2a on the DM distribution obtained as described in appendix C. In the left part of the table we use the 
published unmodulated event rates from [4] , whereas for the right part of the table we adopt the preliminary 
results on surface events contamination at low energies from [21]. 

Appendix C: Alternative procedure for optimizing the bound in presence of un- 
known background 

In sec. 5 we have outlined a method based on a least-square minimization ("procedure 2") 
to find "optimal" values for the coefficients Ui parameterizing the background contribution 
in each bin. Here we provide an alternative method for this task. 

Let us consider the inequalities (38) and (39) as a system of equations for the Ui. We 
observe that they have a triangular structure with respect to uji. the inequality for a given 
index i depends only on ujj with j > i, i.e., the inequality for i = depends only on un, 
ioT i = N — 1 it depends on un-i, i^n, and so on. From the fact that a^, (k)^ , ji are positive 
follows that the bound for a given i is weakest if Ui is as large as possible and uj with j > i 
are as small as possible. We can implement this requirement in an iterative way: replacing 
the inequality sign by an equality, we obtain a system of linear equations in Ui. We can 
solve this system by starting at i = N and going up step by step to i = 1, in each step 
obtaining a value for the corresponding Ui, to be used in the next step. If the solution for 
Ui is less than one, it will be the smallest allowed value, and hence the bound for i — 1 will 
be weakest. If the solution for Ui is larger than one the bound is violated by bin i and we 
have to set the corresponding Ui = 1 (the largest physically allowed value). In that way we 
obtain a set of Ui corresponding to the most conservative choice of background contributions 
to the unmodulated rate. 

This method treats the inequalities as "hard bounds" and does not take into account the 
fact that the amplitudes are subject to experimental uncertainties. A conservative way to 
include uncertainties is to apply this procedure as described above but using instead of the 
l.h.s. of (38) and (39) its central value minus n standard deviations. In this way one aims to 
satisfy the inequalities within the na interval for each bin. We show the results of such an 
analysis for CoGeNT data in tab. III. Typically one finds qualitatively similar but slightly 
stronger bounds than for procedure 1, although if the method is considered at higher CL 
bounds may become weaker or even disappear, as visible in the table for the columns at 
95% CL. Note, however, that subtracting na from the mean value for all bins leads to very 
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conservative limits. 
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